Astronomy & Astrophysics manuscript no. HD10180 


©ESO 2010 


November 24, 2010 





The HARPS search for southern extra-solar planets^ 

XXVII. Up to seven planets orbiting HD 10180: 
probing the architecture of low-mass planetary systems 

C. Lovis\ D. Segransan\ M. Mayor\ S. Udry\ W. Benz^, J.-L. Bertaux^, F. Bouchy'^'^, A. C. M. Correia*', 
J. Laskar'^, G. Lo Curto^ C. Mordasini''-2, F. Pepe^, D. Queloz\ and N. C. Santosi°'i 

' Observatoire de Geneve, Universite de Geneve, 51 ch. des Maillettes, CH-1290 Versoix, Switzerland 

e-mail: christophe . lovisOunige . ch 
^ Physikalisches Institut, Universitat Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland 

^ Universite Versailles Saint-Quentin, LATMOS-IPSL, 11 Boulevard d'Alembert, F-78280 Guyancourt, France 
Institut d' Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie, 98bis Bd Arago, F-75014 Paris, France 

' Observatoire de Haute-Provence, CNRS/OAMP, F-04870 St. Michel 1' Observatoire, France 
Department of Physics, I3N, University of Aveiro, Campus Universitario de Santiago, 3810-193 Aveiro, Portugal 
ASD, IMCCE-CNRS UMR8028, Observatoire de Paris, UPMC, 77 Av. Denfert-Rochereau, F-75014 Paris, France 
European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei Miinchen, Germany 
Max-Planck-Institut fUr Astronomic, Konigstuhl 17, D-69117 Heidelberg, Germany 

Departamento de Fisica e Astronomia, Faculdade de Ciencias, Universidade do Porto, 4150-762 Porto, Portugal 



Received 12 August 2010 / Accepted 18 November 2010 

ABSTRACT 

Context. Low-mass extrasolar planets are presently being discovered at an increased pace by radial velocity and transit surveys, 
opening a new window on planetary systems. 

Aims. We are conducting a high-precision radial velocity survey with the HARPS spectrograph which aims at characterizing the 
population of ice giants and super-Earths around nearby solar-type stars. This will lead to a better understanding of their formation 
and evolution, and yield a global picture of planetary systems from gas giants down to telluric planets. 

Methods. Progress has been possible in this field thanks in particular to the sub-ms"' radial velocity precision achieved by HARPS. 
We present here new high-quality measurements from this instrument. 

Results. We report the discovery of a planetary system comprising at least five Neptune-like planets with minimum masses ranging 
from 12 to 25 M^, orbiting the solar-type star HD 10180 at separations between 0.06 and 1.4 AU. A sixth radial velocity signal is 
present at a longer period, probably due to a 65-Me object. Moreover, another body with a minimum mass as low as 1.4 Me may 
be present at 0.02 AU from the star. This is the most populated exoplanetary system known to date. The planets are in a dense but 
still well-separated configuration, with significant secular interactions. Some of the orbital period ratios are fairly close to integer 
or half-integer values, but the system does not exhibit any mean-motion resonances. General relativity effects and tidal dissipation 
play an important role to stabilize the innermost planet and the system as a whole. Numerical integrations show long-term dynamical 
stability provided true masses are within a factor ~3 from minimum masses. We further note that several low-mass planetary systems 
exhibit a rather "packed" orbital architecture with little or no space left for additional planets. In several cases, semi-major axes are 
fairly regularly spaced on a logarithmic scale, giving rise to approximate Titius-Bode-like (i.e. exponential) laws. These dynamical 
architectures can be interpreted as the signature of formation scenarios where type I migration and interactions between protoplanets 
play a major role. However, it remains challenging to explain the presence of so many Neptunes and super-Earths on non-resonant, 
well-ordered orbits within ~l-2 AU of the central star. Finally, we also confirm the marked dependence of planet formation on both 
metallicity and stellar mass. Very massive systems are all found around metal-rich stars more massive than the Sun, while low-mass 
systems are only found around metal-deficient stars less massive than the Sun. 

Key words. Planetary systems - Stars: individual: HD 10180 - Techniques: radial velocities - Techniques: spectroscopic 



1. Introduction 

Over the past 15 years, the field of extrasolar planets has been 
witnessing uninterrupted developments and several major mile- 
stones. Among these one can mention: the initial proo f of exis- 
tence of extrasolar gas giants dMavor & Ouelodll995h . the first 
detection of a transiting planet, providing a precise mass and 
radius ( Charbonneau et al. 2000; Henrv et al. 2000), the detec- 



* Based on observations made with the HARPS instrument on the 
ESO 3.6-m telescope at La Silla Observatory (Chile), under program 
IDs 072.C-0488 and 183.C-0972. 



tion of a large sample of gas giants with a variety of masses 
and orbital properties, the characterization of bulk properties 
and atmospheres of transiting gas giants, and the detection of 
objects in the Neptune-mass and super-Earth range. Recently, 
two transiting planets with masses and ra dii close to those 
of Earth have been discovered: CoRoT-7b (iLegeret al.l 120091: 
lOueloz etal.1 12009) and GJ 1214 b dCharbonneau e t al. 2009. 
High-precision radial velocity surveys are now able to find plan - 
ets with minimum masses as low as 1 .9 Me (Ma vor et al.l2009ah . 
Preliminary results from the HARPS survey are hinting at a 
large population of Neptune-like objects and super-Earths within 
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~0.5 AU of solar-type stars (iLovis et alJl2009l) . Moreover, hun- 
dreds of small-radius candidate planets have been announced 
by the Kepler Team (Borucki & the Kepler Team 2010). Clearly, 
the exploration of the low-mass planet population has now fully 
started, and will become the main focus of the field in the coming 
years. It is expected that the characterization of planetary sys- 
tem architectures, taking into account all objects from gas giants 
to Earth-like planets, will greatly improve our understanding of 
their formation and evolution. It will also allow us to eventu- 
ally put our Solar System into a broader context and determine 
how typical it is in the vastly diverse world of planetary sys- 
tems. The characterization of a significant sample of low-mass 
objects, through their mean density and some basic atmospheric 
properties, is also at hand and will bring much desired insights 
into their composition and the physical processes at play during 
planet formation. 

As part of this broad effort to explore the low-mass exoplanet 
population, we are conducting a high-precision radial velocity 
survey of about 400 bright FG K stars in the solar n eighbourhood 
using the HARPS instrument dMavor et al.ll200 3V Observations 
of this sample were obtained during HARPS GTO time from 
2004 to 2009 (PI: M. Mayor), and then were continued as an 
ESO Large Program (PI: S. Udry) until today. Several stars from 
this survey h ave already revealed orbiting low-m ass objects: 
HP 160691 (ISantos et al.' '2004; 'Pepe et al."2007'), HD 69830 
2006a), HD 4308 (Udry et al. 2006), HD 40307 
2009bl) . HD 47186, HD 1814 33 (iBouchv et all 



Table 2. Fundamental properties of HD 10180. 



(Lovis et al.1 



dMavor et al 



120091) and HD 90156 (iMordasini et all 1201 Ol) . More and more 
candidates are detected as measurements accumulate, and many 
new systems are about to be published (Queloz et al., Udry et 
al., Segransan et al., Benz et al., Dumusque et al., Pepe et al., in 
prep.). Following 400 stars to search for radial velocity signals 
at the ms ' level requires a lot of telescope time, and this sur- 
vey is by construction a long-term project. Over the years, we 
chose to focus on a smaller sample of stars showing a low level 
of chromospheric activity to minimize the impact of stellar noise 
on our planet detection limits. Based on measured Call activity 
levels log/?[jj., we kept about 300 stars which we are monitor- 
ing regularly. Once a sufficient number of observations has been 
gathered for each star, we will be able to derive important statis- 
tical properties of the low-mass planet population (Mayor et al., 
in prep.). 

In this paper we present the discovery of a new low-mass 
planetary system comprising at least 5 Neptune-mass planets 
and, probably, a longer-period object and a close-in Earth-mass 
planet. The parent star is the GIV dwarf HD 10180, located 39 
pc away from the Sun towards the southern constellation Hydrus. 



2. Observations and data reduction 

The data presented in this paper have been obtained with the 
HARPS spectrograph at the ESO 3.6-m telescope at La Silla 
Observatory (Chile). This instrument has demonstrated a sub- 
ms~' radial veloc i ty precision over m ore than 6 years (e.g. 
iLovis et al.l[2006bt iMavoret al.ll2009bh and has led to the de- 
tection of the majority of the low-mass planets known to date. 

We have obtained a total of 190 data points on HD 10180, 
spread over more than 6 years. This star is part of the high- 
precision planet-search sample of about 400 stars that we have 
been following closely since 2004. Exposure times were set to 15 
min to average out stellar oscillations. The achieved SNR at 550 
nm ranges from 120 to 270, depending on weather conditions. 
The estimated photon noise limit to the radial velocity precision 
ranges from 80 to 30 cms"', respectively. Including other mea- 



Parameter 




HD 10180 


Spectral type 




GIV 


V 


[mag] 


7.33 


B-V 


[mag] 


0.629 


n 


[mas] 


25.39 ± 0.62 


Mv 


[mag] 


4.35 




[K] 


5911 ± 19 


logg 


[cgs] 


4.39 ± 0.03 


[Fe/H] 


[dex] 


0.08 ± 0.01 


L 


[Le] 


1.49 ± 0.02 


M. 


[Mo] 


1.06 ± 0.05 


V sin I 


[kms-'] 


<3 






-5.00 




[days] 


24 ± 3 


Age (logi?;^,^) 


[Gyr] 


4.3 ± 0.5 



surable instrumental errors (wavelength calibration, noise on in- 
strumental drift measurement), we obtain error bars between 1.3 
and 0.4ms '. This does not include other instrumental system- 
atics like telescope guiding (light injection) errors, which are ex- 
pected to be small but difficult to estimate. Data reduction was 
performed with the latest version of the HARPS pipeline (see 
Lovis et al. 201 1, in prep, for a more detailed description). 

The end products of the reduction are barycentric radial ve- 
locities with internal error bars, bisector span measurements, 
parameters of the cross-correlation functions (FWHM and con- 
trast), and Call activity indices S and logR'^^^. Being a stabi- 
lized, well-controlled instrument, HARPS provides not only pre- 
cise radial velocities, but also precise spectroscopic indicators in 
general, which is very useful to better understand the stars un- 
der consideration ( see e.g. the case of the active star CoRoT-7, 
lOueloz et al.ll2009l) . The whole set of radial velocities and spec- 
troscopic measurements of HD 10180 can be found in Table [T] 
and in electronic form at CDS. 



3. Stellar properties 

The fundamental properties of HD 10180 (GIV, y=7.33) are 
taken from the Hi pparcos catal ogue (ESA 1997) and the spectro- 
scopic analysis bv lSousa et al.l (120081) . They are listed in Table|2] 
HD 10180 is a solar-type star with a mass M = 1.06 ± 0.05 Mq 
and metallicity [Fe/H] - 0.08 + 0.01 dex. With a mean activity 
index log/?[jj, - -5.00, measured on the HARPS spectra pre- 
sented here, it is clearly an inactive star. Furthermore, it does not 
show any well-defined activity cycle such as the solar one (the 
rms dispersion of the logR'^^^ measurements is only 0.012 dex). 
Given this low activity level and the early-G spectral type, we 
expect a stellar RV jitter at the level of ~1 ms ' based on com- 
parisons with similar stars in the HARPS sample (see Dumusque 
et al. 2011, in prep., Lovis et al. 2011, in prep.). Among inac- 
tive stars, early-G dwarfs appear to have slightly more jitter than 
late-G and early-K dwarfs, possibly due to more vigorous con- 
vection and thus increased granulation noise. We thus adopt a 
value of 1.0 ms"' for the stellar jitter in this paper and add this 
number quadratically to the instrumental error bars. The main 
purpose for doing this is to avoid large, unjustified differences in 
the individual weights (vv,= l/crf) used in the -minimization 
process. 
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Fig. 1. Close-up views of the window function of the radial ve- 
locity measurements, centered on the regions that have signifi- 
cant peaks and that may induce aliases in the RV data. 



4. Analysis of the radial velocity data 

4.1. Detection of 5 strong signals 

The raw rms dispersion of the radial velocities is 6.42 m s ' , well 
above instrumental errors and the expected stellar jitter, hint- 
ing at the presence of planets orbiting HD 10180. We proceed 
to an analysis of the data using the Yorbit software, an analy- 
sis package for radial velocity, astrometric and photometric data 
featuring a genetic algorithm and various tools for exoplanet 
search (Segransan et al. 2011, in prep.). We perform succes- 
sive multi-Keplerian fits to the RV data, and use Lomb-Scargle 
periodograms to search for periodic signals of significant am- 
plitude in the residuals to the successive models. We use the 
generalized Lomb-Scargle (G LS) periodogram as described in 
IZechmeister & Kurstei^('2009). False-alarm probabilities (FAPs) 
are computed by performing random permutations of the data, 
keeping the observing times fixed, and computing their peri- 
odograms. For each trial the peak power is recorded, and then 
the power of the real signal is compared to the peak power distri- 
bution of the permuted datasets. Each time a signal is considered 
significant, we include it in the multi-Keplerian model and pro- 
ceed further, assuming the radial velocity signals are due to or- 
biting planets. To do so, a genetic algorithm is run to efficiently 
explore the parameter space around suspected orbital periods. 
Once the population of solutions has converged towards the sin- 
gle best-fitting region of parameter space, a final Levenberg- 
Marquardt minimization is performed to reach the deepest 
minimum found. 

Identification of signals in periodograms is sometimes am- 
biguous due to the presence of aliases, induced by the non- 
random time sampling of the observations. This problem has 
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Fig. 2. Successive GLS periodograms of the HD 10180 RV time 
series, where the main signal (indicated by an arrow) is sub- 
tracted at each step by adding one more planet in the multi- 
Keplerian model. FAP thresholds of 10%, 1% and 0.1% are in- 
dicated as dashed lines. 
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Fig. 3. Peak power distributions of the random permutations of the residuals to the successive multi-Keplerian models. The peak 
power of the actual data is shown as a full vertical line, while 10%, 1% and 0.1% thresholds are shown as dashed lines. From left 
to right and top to bottom, the peak power corresponds to periods of 50 d, 5.8 d, 122 d, 2200 d, 16 d, 600 d and 1.18 d. The first 5 
signals have extremely high significance, while the 6th and 7th signals have 0.6% and 1.4% FAP, respectively. Nothing significant 
remains after subtracting the 7-Keplerian model. 



been long known in several fields of astrophysics dealing with 
ground-based time series analysis, and is also relevant to radial 
velocity searches for exoplanets (see Dawson & Fabrvckv 2010, 
for a recent overview). Aliases occur at frequencies separated 
from the true peak by a frequency difference at which the win- 
dow function shows significant power. Fig. [T] shows 3 close-up 
views of the window function of the HD 10180 measurements, 
centered around f - Q, f - Id"' and / = 2 d"'. These are the 
regions which contain significant peaks and which could con- 
tribute aliases in the frequency range of interest. As expected, the 
most prominent features are found at / = 0.0027 d"' (1 year), / 
= 1.0027 d-i (1 sidereal day), /= 1.0000 d-i (1 day), / = 2.0027 
d ' and / = 2.0055 d In the following we will pay particular 
attention to possible aliases induced by these peaks. 

Fig.|2]shows the successive GLS periodograms of the radial 
velocity data. At each step, the main peak is identified and con- 
sidered significant if its false-alarm probability (FAP) is smaller 
than 10"^. FAP thresholds of 10 ', 10"^ and 10 ai-e indicated 
by horizontal lines in Fig. |2l The peak power distributions of 
the shuffled datasets are shown in Fig. [3] As can be seen from 
both figures, 5 very significant signals can be successively found 
in the data, at periods 50 d, 5.8 d, 122 d, 2200 d and 16 d. All 
the corresponding peaks are already clearly apparent in the pe- 
riodogram of the raw data. Each of them has a FAP far below 
10"^, as can be seen from the peak power distributions. We suc- 
cessively subtract these signals by performing multi-Keplerian 
fits to the data, whose characteristics can be found in Table |3] 
The final, 5-Keplerian fit to the data yields periods of 5.760 d, 
16.35 d, 49.74 d, 122.4 d and 2231 d for the 5 signals. 

Various aliases of these signals are present in the peri- 
odograms. It turns out that the 5 successive highest peaks are 
always located at periods larger than 2 d, and in each case two 
other peaks potentially corresponding to their 1-d aliases are 
clearly seen. Given a highest-peak frequency /o, these are found 



at |/o + 1.0027 d I and show amplitudes similar to each other 
but lower than the main peak. Given the properties of the win- 
dow function, this is the expected pattern if the true signal is the 
lower-frequency peak at P > 2 d. On the contrary, if the true 
signal was one of the two peaks near 1 d, the peak at P > 2 d 
would be a 1 -d alias and the other high-frequency peak would 
be a 0.5-d alias (2.0055 d"'). In this case, one would expect a 
large difference in amplitude between the two peaks near 1 d, 
and the low-frequency peak should be of intermediate strength. 
We checked the global pattern of aliases on the periodogram of 
the raw data for the 5 strong signals, extending the computation 
to frequencies around 2d'. For each signal, we verified that 
the amplitude envelope outlined by the relevant peaks (the main 
peak plus its supposed 1-d and 0.5-d aliases) corresponds to the 
amplitude envelope in the window function. In each case we find 
a symmetric amplitude pattern centered on the lower-frequency 
peak and thus we conclude that, for the 5 strong signals seen 
in the raw data, the correct periods are the "long" ones, and the 
forest of peaks around P = 1 d are aliases. 

In summary, we obtain a 5-Keplerian fit with periods 5.760 d, 
16.35 d, 49.74 d, 122.4 d and 2231 d. Table[3]Hsts the statistical 
characteristics of the successive multi-Keplerian models. 

4.2. The 600-d signal 

After subtraction of these 5 signals, the periodogram of the resid- 
uals still contains appreciable power, with peaks at P = 600 
d, 227 d and around 1 d. One immediately sees that the first 
two peaks are aliases of each other with the 1-year frequency 
(1/600 + 1 /365 ^ 1 /227). The FAP of the highest peak (P = 600 
d) is 0.6%, and we thus consider it as significant. Given that the 
spectral window has relatively high sidelobes at the 1 -year fre- 
quency, it is not surprising that a fraction of the power is leaking 
into the 1-year aliases. Here we assume that the correct period is 
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Table 3. Characteristics of the successive multi-Keplerian mod- 
els fitted to the data. P is the period of the last signal added to 
the model. 



Model 


P 




T 


T 

xi 


cr(O-C) 




[days] 








[m s-'] 


constant 




1 


5741.6 


30.38 


6.42 


kl 


50 


6 


4065.9 


22.10 


5.45 


k2 


5.8 


11 


2568.5 


14.35 


4.29 


k3 


122 


16 


1540.1 


8.85 


3.27 


k4 


2200 


21 


824.0 


4.88 


2.36 


k5 


16 


26 


356.4 


2.17 


1.57 


k6 (£5=0) 


600 


29 


276.0 


1.71 


1.39 


k7(ei=e6=0) 


1.18 


32 


237.1 


1.50 


1.27 



not one of the peaks around 1 d because it is highly unlikely to 
find a planet with an orbital period so close to 1 d (closer than 
0.005 d). To determine which period (600 d or 227 d) is the most 
likely one, we proceed in the following way. We perform simu- 
lations in which we take the residuals from the full 6-Keplerian 
model (with either P - 600 d or P = 227 d as the 6th signal), 
shuffle them using permutations, inject an artificial signal at ei- 
ther P - 600 d or P = 227 d, and compare the overall pattern 
of peaks in the resulting periodograms to the observed one. The 
periods, amplitudes and phases of the injected signals are drawn 
from Gaussian distributions centered on the fitted values in the 
6-Keplerian model. As a quantitative measure of the similarity 
of periodograms, we check whether the 3 highest peaks are the 
same as in the actual data. As actual data to compare to, we take 
the residuals to the 6-Keplerian model with the fitted 6th signal 
added (either at P = 600 d or P = 227 d). In this way we compare 
the simulated data to the 5-Keplerian residuals that are closest to 
reality under the two assumptions (the 600-d or the 227-d peak 
is the correct one). We take all these precautions because we 
are dealing with data on which 5 signals have already been sub- 
tracted, and the exact choice of parameters for these 5 signals 
has a significant impact on the alias pattern of the 6th signal. 

The results of the simulations are as follows: when injecting 
a 600-d signal, the 3 highest peaks seen in the actual data (600 
d, 0.9956 d and 0.9989 d) are correctly reproduced in 53% of 
the simulated periodograms. When injecting a 227-d signal, the 
3 highest peaks (227 d, 1.0017 d and 0.9956 d) are recovered in 
only 1.3% of the cases. We thus conclude that a 600-d signal is 
much more likely to correctly reproduce the data than a 227-d 
signal. 

We proceed to fit a 6-Keplerian model to the data, which 
yields;^'^ = 1.71 and cr(O-C) - 1.39 ms"', compared to =2.17 
and o"(0-C) = 1.57ms ' for the 5-Keplerian model. We check 
whether the eccentricity for the 6th, lowest-amplitude signal is 
constrained by the data by fixing it to zero, refitting and looking 
at the reduced;\f^. With;if^ = 1.71, as before, the addition of more 
free parameters in the model is not justified and we adopt the 
zero-eccentricity solution. Finally, it is also worth mentioning 
that the 6-Keplerian solution with P = 600 d has a slightly better 

(1-71) than the corresponding one with P = 227 d (1.82), 
reinforcing the case for the longer period. 

4.3. A potential 7th signal 

Continuing the process one step further, we notice two more 
fairly strong peaks in the periodogram of the residuals to the 
6-Keplerian model (see Fig. |2]l. These are located at periods P 
= 6.51 d and P = 1.178 d. Again, one of them is clearly the 
alias of the other one, this time with the 1 sidereal day period 



(11/6.51 - 1.00271 ^ 1/1.178). After 10,000 random permuta- 
tions of the residuals, we obtain a FAP of 1.4% for the higher 
peak (P = 6.51 d). This is slightly too high to confidently claim 
a detection, but it is nevertheless intriguing. We proceed to fit 
this possible 7th signal, but first the correct period has to be de- 
termined. As before, we perform simulations by injecting artifi- 
cial signals into the 7-Keplerian residuals at both periods. When 
signals are added at P = 1.178 d, the 3 highest peaks found in 
the actual data (6.51 d, 1.178 d, 1.182 d) are reproduced in only 
0.2% of the cases. With signals injected at P = 6.51 d, the 3 high- 
est peaks (6.51 d, 1.178 d and 1.182 d, as before) are recovered 
in just 0.7% of the simulated periodograms. This slightly favors 
the 6.5 1-d period, but only marginally so. Above all, the very 
low "success" rate of the simulations seems to indicate that this 
method cannot reliably decide which peak is the correct one. We 
also note that as much as 1 1 % of the simulated periodograms 
with an injected signal at P = 1 . 178 d yield a highest peak at P = 
6.51 d. It would therefore be very speculative to draw any con- 
clusion from the fact that the observed highest peak is at P = 6.5 1 
d. In summary, it is likely that the present data are simply not suf- 
ficient to distinguish between aliases. The strong 1-d aliases in 
the window function (86% of the main peak) are the main obsta- 
cles to overcome, which means that several data points spread 
within nights will be required in the future to resolve this issue. 

In the meantime, we use another kind of argument: if this 
7th signal is caused by an orbiting planet, then it is very unlikely 
that the system would be dynamically stable with two objects 
at P = 5.76 d and P = 6.51 d, especially considering the 13- 
Me minimum mass of the former planet (see Sect. |6]for further 
discussion on this point). We thus adopt P = 1.178 d as the only 
viable period from a dynamical point of view. The 7-Keplerian 
model, with eccentricities of the two lowest-amplitudes signals 
fixed to zero, has;^'^ = 1.50 and cr(O-C) = 1.27 ms"'. 

As a last step, an inspection of the periodogram of the resid- 
uals to the 7-Keplerian model does not reveal any peak with a 
FAP below 10%, thus ending the search for signals in the RV 
data. 

In summary, we firmly detect 6 signals with periods 5.760 
d, 16.36 d, 49.75 d, 122.7 d, 602 d and 2248 d. A 7th signal 
with P = 1.178 d has FAP 1.4% and requires confirmation. The 
6-Keplerian model has 29 free parameters for 190 data points, 
andx^ = 1.71. This value goes down to 1.50 for the 7-Keplerian 
model. The rms dispersion of the residuals is 1.39ms"', down 
to 1.27 ms"' in the 7-Keplerian model. We thus have a good 
fit to the data, confirming that the adopted stellar jitter value 
(1.0ms"') is reasonable. 

5. Origin of the radial velocity signals 

So far we have assumed that all RV signals in the data are caused 
by orbiting planets. Obviously, with the small semi-amplitudes 
in play (between 0.81 and 4.5 ms"'), this assumption should be 
further verified. For this purpose we study the behaviour of sev- 
eral spectroscopic indicators: bisector velocity span, FWHM of 
the cross-correlation function and Call activity index logP^j^. 
Time series of these indicators and the corresponding GLS peri- 
odograms are shown in Fig.|4] 

5.1. Bisector velocity span and CCF FWHM 

The global rms dispersion of the bisector span is only 1.33ms"', 
a remarkable stability over more than 6 years. An analysis of 
the bisector periodogram reveals no significant power at any fre- 
quency, indicating a very quiet star In particular, no power is 
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Fig. 5. Long-period radial velocity signal (P - 2248 d) as ob- 
tained from the 7-Keplerian model, as a function of logR'^^. A 
weak correlation seems to be present between the two quantities 
(R = 0.44). 



seen around the estimated rotation period (~24 d) or at any of 
the 7 frequencies seen in the radial velocities. The same is true 



for the CCF FWHM, although some more structure seems to 
be present in the data. The FWHM periodogram shows no peak 
withaFAP below 3%. 



5.2. Call activity index 

We now turn to the activity index logR'^^^. We see that a long- 
period modulation is present in these data, although of very low 
amplitude (a few 0.01 dex) and at the same level as the short- 
term scatter This behaviour is not typical of Sun-like magnetic 
cycles, which have long-term amplitudes an order of magnitude 
larger. The peak-to-peak variations here are only 0.06 dex, com- 
pared to 0.2-0.3 dex for typical activity cycles in solar-type stars. 
Also, the data seem to show a modulation only in the second 
half of the observing period. Thus, HD 10180 does not presently 
show an activity cycle like the Sun and, with an average log R'^^^ 
of -5.00 over 6 years, seems to be in a very quiet state. Its fun- 
damental parameters do not indicate a subgiant status, which 
could have explained the low and stable activity level. More 
likely, HD 10180 either has a magnetic cycle with a period much 
longer than 10 years, or it is in a relatively quiet phase of its 
main-sequence lifetime, with sporadic, weak variations in activ- 
ity level. If the latter explanation is true, we might be witnessing 
an activity state similar to the Maunder minimum of the Sun in 
the XVIIth century. 

As shown by Lovis et al. (2011, in prep.), magnetic cycles 
do induce RV variations in solar-type stars, at a level that de- 
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pends on spectral type. Studying a large sample of stars observed 
with HARPS, they were able to measure the degree of correla- 
tion between activity level and radial velocities. It turns out that 
the sensitivity of RVs to activity variations increases with in- 
creasing temperature, early-G dwarfs being much more affected 
than K dwarfs, wh ich are almost immune to this phenomenon 
(ISantos et al.l l2010). In the case of HD 10180, the iogR'.^^ peri- 
odogram exhibits a very significant peak around 2500 d, as can 
already be guessed by eye from the time series (see Fig. |4]i. A 
comparison with the raw radial velocity curve shows that the 
long-period RV signal at P ~ 2200 d might be correlated to the 
log7?[jj^ signal. Fig. |5] shows this long-period RV signal (plus 
residuals), as obtained from the 7-Keplerian model, as a func- 
tion of logR'^^. The Pearson's weighted correlation coefficient 
R is 0.44, indicating a weak positive correlation between the 
two quantities. The measured slope of the RV-log/?^j^ relation 
is 0.92 ± 0.13ms"' per 0.01 dex. This number, although not 
well constrained, is compatible with the predicted sensitivity to 
activity of HD 10180. Indeed, for T^g = 5911 K, the activity- 
RV relation of Lovis et al. (2011, in prep.) gives 0.99 ms"' per 
0.01 dex. It is thus possible that the long-period RV signal is not 
of planetary origin, but the result of the varying fraction of the 
stellar surface covered by magnetic regions over time. However, 
HD 10180 does not show a typical solar-type magnetic cycle, 
and the quality of the correlation is quite poor For example, the 
activity data around JD=53700 exhibit a large scatter instead of 
a uniformly low value, as would be expected from the RV val- 
ues in the case of an activity-induced signal. Moreover, the CCF 
FWHM does not show a clear correlation with log R{i^, in con- 
trast to the stars studied by Lovis et al. (201 1, in prep.). Also, the 
fitted semi-amplitude of the RV signal is 3.11ms"' while the 
fitted semi-amplitude on the logR'^^ data is 0.011 dex, which 
appears to be a factor ~3 too low to account for the RV signal 
considering a sensitivity of ~1 ms"' per 0.01 dex. Such a fac- 
tor is large compared to the scatter in the Lovis et al. (2011, in 
prep.) relation and it thus seems unlikely that HD 10180 could 
be so sensitive to activity. 

In conclusion, we favor the planetary interpretation for this 
long-period RV signal, although some doubt remains on its ori- 
gin. Future observations will likely solve this issue. In particu- 
lar, it will be interesting to see whether radial velocities follow a 
downward trend in the near future, as would be expected in the 
planetary case, or whether they follow the more chaotic varia- 
tions of the activity index. 

There is no indication that the 6 other RV signals might be 
due to stellar activity. None of them are related to the expected 
stellar rotation period (24 ± 3 d) or to its harmonics. We tried 
to detect the rotation period in the log R[^^ data, but found no 
convincing signal. This confirms that HD 10180 is a very quiet 
star. Another strong argument in favor of the planetary interpre- 
tation is the very high significance of the signals in the peri- 
odograms. It is a clear indication that the signals are coherent 
in time, as expected from the clock-like signature of an orbiting 
planet. Activity-related phenomena like spots and plages have 
short coherence times in inactive stars (of the order of the ro- 
tation period), and are not able to produce such well-defined, 
high-significance peaks in a dataset spanning more than 6 years. 
We thus conclude that a planetary origin is the only viable inter- 
pretation for the 6 RV signals between 1 d and 600 d. 

6. The HD 10180 planetary system 

From now on we will assume that all 7 RV signals are real and 
of planetary origin. We thus keep the 7-Keplerian model - and 
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Fig. 6. Radial velocity time series with the 7-Keplerian model 
overplotted. The lower panel shows the residuals to the model. 



the fit - as they were obtained in Sect.|4] Table|4]gives the orbital 
parameters of all planets in the system. Fig.|6]shows the full ra- 
dial velocity curve as a function of time, while the phased RV 
curves for all signals are shown in Fig. |7] Uncertainties on the 
fitted parameters have been obtained with Monte Carlo simula- 
tions, where actual data points are modified by drawing from a 
Gaussian distribution with a standard deviation equal to the error 
bar on each point. The modified datasets are then re-fitted, ini- 
tializing the non-linear minimization with the nominal solution, 
and confidence intervals are derived from the obtained distribu- 
tions of orbital elements. 

The orbital solution given in Table |4] can be further im- 
proved by using N-body integrations and adding constraints on 
the orbital elements obtained from dynamical considerations 
(see Sect. |7]l. To decouple things, we first show here the solu- 
tion obtained with a multi-Keplerian model and based on the RV 
data alone. Its main advantage is its much higher computational 
simplicity, which allows us to derive reliable error bars through 
Monte Carlo simulations. 

The HD 10180 planetary system is unique in several re- 
spects. First of all, it presents five Neptune-like planets, orbiting 
between 0.06 and 1.4 AU from the central star With minimum 
masses between 12 and 25 M®, this represents a relatively large 
total planetary mass within the inner region of the system, and 
multi-body migration processes are likely needed to explain this. 
Besides these Neptune-mass objects, the system has no massive 
gas giant. At most, it contains a small Saturn (msin/ = 65 M®) 
at 3.4 AU. In fact, the present data show no detectable long-term 
drift and allow us to exclude any Jupiter-mass planet within ~10 
AU for an edge-on system. At the other extreme of the mass 
and semi-major axis scales, the system also probably contains 
an Earth-mass object (msin / =1.4 M®) orbiting at only 0.022 
AU. This is the planet with the lowest minimum mass found to 
date, and may represent another member of a hot rocky planet 
population that is starting to emerge (e.g. CoRoT-7b, GJ 581 e, 
HD 40307 b). The discovery of the HD 10180 planets highlights 
once again how diverse the outcome of planet formation can be. 
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Table 4. Orbital and physical parameters of the planets orbiting HD 101 80, as obtained from a 7-Keplerian fit to the data. Error bars 
are derived using Monte Carlo simulations. A is the mean longitude (A-M + a))at the given epoch. 



Parameter 


[unit] 


HD 10180 b 


HD 10180 c 


HD 10180 d HD 10180 e HD 10180 f 


HD 10180 g 


HD 10180 h 


Epoch 
i 

V 


[BJD] 
[deg] 
[km s-'] 






2,454,000.0 (fixed) 

90 (fixed) 
35.53015 (±0.00045) 








P 
A 
e 

0) 

K 


[days] 
[deg] 

[deg] 

[ms^'j 


1.17765 
(±0.00018) 
186 
(±11) 
0.0 
(fixed) 

0.0 
(fixed) 
0.81 
(±0.15) 


5.75962 
(±0.00028) 
238.1 
(±1.9) 
0.077 
(±0.033) 
279 

/+63x 
'•-25'' 

4.54 
(±0.15) 


16.3567 
(±0.0043) 

197.9 
(±3.3) 

0.143 
(±0.058) 

292 

/+42x 
"•-10-' 

2.93 
(±0.16) 


49.747 
(±0.024) 
102.3 
(±2.2) 
0.065 
(±0.035) 
174 

'•-59'' 

4.25 
(±0.18) 


122.72 
(±0.20) 
251.4 
(±3.2) 
0.133 
(±0.066) 
265 

/•-f78 ^ 
'•-209^ 

2.95 
(±0.18) 


602 
(±11) 

321 
(±11) 

0.0 
(fixed) 

0.0 
(fixed) 

1.56 
(±0.21) 


2248 

'-106'' 

238.8 
(±4.1) 
0.151 
(±0.072) 
194 

3.11 
(±0.21) 


m sin ; 
a 


[Mel 
[AU] 


1.39 
(±0.26) 
0.02225 
(±0.00038) 


13.17 
(±0.63) 
0.0641 
(±0.0011) 


11.94 
(±0.75) 
0.1285 
(±0.0022) 


25.4 
(±1.4) 
0.2699 
(±0.0045) 


23.6 
(±1.7) 
0.4928 
(±0.0083) 


21.4 
(±3.0) 
1.423 
(±0.030) 


65.3 
(±4.6) 
3.42 

^-0 n> 


N 

Span 
rms 

x; 


[days] 
[ms-'] 








190 
2428 
1.27 
1.50 









7. Dynamical analysis 



With such a complex system at hand, various dynamical stud- 
ies are in order. For the first time, an extrasolar planetary sys- 
tem comes close to the Solar System as far as the number of 
bodies involved is concerned. An obvious aspect to be checked 
is whether the fitted orbital solution is dynamically stable over 
Gyr timescales (the age of the star). The stability of such a sys- 
tem is not straightforward, in particular taking into account that 
the minimum masses of the planets are of the same order as 
Neptune's mass and the fitted eccentricities are relatively high 
when compared with the eccentricities of the planets in the Solar 
System. As a consequence, mutual gravitational interactions be- 
tween planets in the HD 10180 system cannot be neglected and 
may give rise to instability. This said, it must be noted that the 
fitted eccentricities of all planets are different from zero by less 
than 2.5 cr (their probability distributions are close to Gaussians), 
thus making it difficult to discuss this issue based on the RV data 
only. Considering the low radial velocity amplitudes induced by 
these objects, it will be challenging to better constrain these ec- 
centricities in the future with the RV method. 



7.1. The secular planetary equations 



Over long times, and in absence of strong mean motion res- 
onances, the variations of the planetary elliptical elements are 
well described by the secular equations, that is the equations ob- 
tain ed after averag ing over the longitudinal motion of the planets 
(see iLaskall 1 9901) . The secular system can even be limited to its 
first order and linear part, which is usually called the Laplace- 



Lagrange system (see lLaskar|[T990l) which can be written using 
the classical complex notation Zk - ejte™' for k - b, c, h 



d_ 
dt 



'Zb- 








= i(L) 




. Zh , 




. Zh , 



(1) 



This linear equation is classically solved by diagonalizing 
the real matrix (L) through the hnear transformation to the 
proper modes 



(2) 



Using the initial conditions of the fit in Table|7] we have com- 
puted analytically the (real) Laplace-Lagrange matrix (L), and 
derived the (real) matrix 5 5 of its eigenvectors which gives the 
relation from the original eccentricity variables to the proper 
modes mj. The matrix 5 5 is given in Table|5] After the transfor- 
mation (l2]i, the linear system ([TJ becomes diagonal 



'Zb^ 




' Ml 




^SS 




. Zh , 




. M7 . 



—(uk) = Diagigi, . 
dt 



(3) 



where gi, . . .,gj are the eigenvalues of the Laplace-Lagrange 
matrix (L). The solution is then trivially given for all A: = 1 , . . . , 7 
by 

= M,(0)e's" , (4) 

while the secular solution is obtained through (|2|. It can be noted 
that as the matrix (L), and thus S S , only depend on the masses 
and semi-major axis of the planets, they do not change much for 
different fits to the data because the periods and masses are well 
constrained (for a given inclination of the system). 

The secular frequencies gk that are responsible for the pre- 
cession of the orbits and for most of the eccentricity variations 
are given in Table|6] 
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0.993728744 
-0.105551000 

0.036827945 
-0.002421676 

0.000090033 
-0.000000038 
-0.000000002 



-0.095494900 
-0.630354828 

0.763187708 
-0.104900945 

0.008592412 
-0.000001984 

0.000000033 



0.054659690 
0.710320992 
0.531792964 
-0.438298758 
0.132463357 
-0.000555340 
-0.000012877 



-0.019554658 
-0.288029266 
-0.347823478 
-0.741307569 
0.496102923 
-0.004533518 
-0.000099417 



0.003718962 
0.063179205 
0.111312700 
0.497011613 
0.857137172 
-0.043581594 
0.000096168 



-0.000102684 
-0.001823646 
-0.003535982 
-0.017867723 
-0.039224786 
-0.986552482 
0.157608762 



0.000014775 
0.000265194 
0.000525407 
0.002723007 
0.006228615 
0.157461056 
0.987501625 



Table 6. Fundamental secular period and frequencies. 



k 


Period 


8k (num) 


8k (ana) 


8k 






yr 


deg yr"' 


deg yr"' 


yr"' 




1 


1029.34 


0.349739 


0.358991 


-3.68x10- 


-07 


2 


1453.39 


0.247696 


0.245229 


-3.40x10- 


-09 


3 


3020.08 


0.119202 


0.118471 


-1.11 xlO- 


-09 


4 


4339.70 


0.082955 


0.079644 


-1.43 xlO- 


-10 


5 


13509.96 


0.026647 


0.025290 


-5.16 xlO 


-12 


6 


61517.43 


0.005852 


0.005581 


-3.93 xlO- 


-15 


7 


473061.76 


0.000761 


0.000663 


-8.14 xlO- 


-17 



Note. The period and gj. (num) are computed numerically from the in- 
tegrated solution of Table[7]through frequency analysis, gt (ana) are the 
corresponding frequencies computed from the Laplace-Lagrange linear 
approximation (Eq[2l(, and gk is the value of the damping term in the 
proper mode amplitudes resulting from the tidal dissipation of planet b 
{Uk = i<t(0)e-«'e«')- 



7.2. Stability of tlie sliort-period planet b 

Despite the proximity of planet b to the star, it undergoes strong 
gravitational perturbations from the remaining bodies, due to 
secular interactions. Even for a model where the initial eccen- 
tricity of the planet b is initially set at zero (TablelH, its eccen- 
tricity shows a rapid increase, which can reach up to 0.8 in only 
1 Myr (Fig. [8^). The inclusion of general relativity in the model 
calms down the eccentricity evolution, but still did not prevent 
the eccentricity of planet b to attain values above 0.4, and planet 
c to reach values around 0.3, which then destabilizes the whole 
system (Fig.lSJ)). 

Most of the variations observed in Fig. [8] are well described 
by the secular system ( I1I2I3I ). In particular, the effect of general 
relativity will be to largely increase the first diagonal terms in 
the Laplace-Lagrange matrix (L), which will increase in zi the 
contribution of the proper mode ui and thus decrease the long 
time oscillations due to the contribution of the modes. 

At this stage, one may question the existence of the inner- 
most short-period planet. However, since the mass of this planet 
should be in the Earth-mass regime, we may assume that it is 
mainly a rocky planet. As a consequence, due to the proximity 
of the star, this planet will undergo intense tidal dissipation, that 
may continuously damp its orbital eccentricity. 

7.3. Tidal contributions 

Using a simplified model (ICorrei al l2009l) the tidal variation of 
the eccentricity is 

e ^ -Knfde)(l - e^)e , (5) 



where Me) = (1 -H 45e^/U + Se* + 685eV224 + 255eV448 + 
25e'"/1792)(l - e^)-^^^^/(l + 3e^ + 3e*/8), and 

21 /M^\lR\^ k2 

Mi, is the mass of the star, m and R are the mass and the radius 
of the planet respectively, k2 is the second Love number, and Q 
the dissipation factor 

As for the Laplace-Lagrange linear system, we can linearize 
the tidal contribution from expression (|5]l, and we obtain for each 
planet k a contribution 

ek = -ykSk with jk = Kktik , (7) 

that is, an additional real contribution on each diagonal term of 
the Laplace-Lagrange matrix i (L) 

Zk = -JkZk , (8) 

which thus adds an imaginary part iy^ to the diagonal terms of 
(L). In fact, since apart from planet b, all planetary masses are 
relatively large (Table|4|i, the dissipation in these planets may be 
small, and we will uniquely consider the tidal dissipation in the 
innermost planet b. It should nevertheless be stressed that all the 
eigenvalues of the matrix (L) will be modified, and will present 
an imaginary part that will induce an exponential term in the 
amplitude of the proper modes (Table |6]l. 

Adopting a value similar to Mars k2/Q - 0.0015 (which is a 
minimal estimation for the dissipation), we can see in Fig.|9]that 
the amplitude of the proper mode mi will be rapidly damped in 
a few tens of Myr, whatever its original value. Over more than 
1 Gyr, the amplitude of the second proper mode U2 is also most 
probably damped to a small value. 

In order to obtain a realistic solution that is the result of the 
tidal evolution of the system, it is thus not sufficient to impose 
that the innermost planets have small eccentricity, as this may 
only be realized at the origin of time (Fig.[8]l. It is also necessary 
that the amplitude of the first proper modes, and particularly ui 
are set to small values. This will be the way to obtain a solution 
with moderate variations of the eccentricities on the innermost 
planets, which will then increase its stability. 

7.4. Orbital solution with dissipative constraint 

As the proper modes of the two innermost planets, ui and U2 are 
damped after about 1 Gyr, we expect them to be small for the 
present observations (the age of the star is estimated to be 4 Gyr, 
Table|2]l. The initial conditions for the HD 10180 planetary sys- 
tem should then take into account the tidal damping. We have 
thus chosen to modify our fitting procedure in order to include a 
constraint for the tidal damping of the proper modes mi and U2. 
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Table 7. Orbital parameters for the seven bodies orbiting HD 10180, obtained with a 8-body Newtonian fit to observational data, 
including the effect of tidal dissipation. Uncertainties are estimated from the covariance matrix, and A is the mean longitude at the 
given epoch (A - M + to). The orbits are assumed co-planar. 



Parameter 


[unit] 


HD 10180 b 


HD 10180 c 


HD 10180 d HD 10180 e HD 10180 f 


HD 10180 g 


HD 10180 h 


Epoch 


[BJD] 








2,454,000.0 (fixed) 










[deg] 








90 (fixed) 








V 


[kms-'] 






35.52981 (±0.00012) 








P 


[days] 


1.17768 


5.75979 


16.3579 


49.745 


122.76 


601.2 


2222 






(±0.00010) 


(±0.00062) 


(±0.0038) 


(±0.022) 


(±0.17) 


(±8.1) 


(±91) 


A 


[deg] 


188 


238.5 


196.6 


102.4 


251.2 


321.5 


235.7 






(±13) 


(±2.3) 


(±3.8) 


(±2.4) 


(±3.6) 


(±9.9) 


(±6.0) 


e 




0.0000 


0.045 


0.088 


0.026 


0.135 


0.19 


0.080 






(±0.0025) 


(±0.026) 


(±0.041) 


(±0.036) 


(±0.046) 


(±0.14) 


(±0.070) 


0) 


[deg] 


39 


332 


315 


166 


332 


347 


174 






(±78) 


(±43) 


(±33) 


(±110) 


(±20) 


(±49) 


(±74) 


K 


[ms-'] 


0.78 


4.50 


2.86 


4.19 


2.98 


1.59 


3.04 






(±0.13) 


(±0.12) 


(±0.13) 


(±0.14) 


(±0.15) 


(±0.25) 


(±0.19) 


m sin ; 


[Me] 


1.35 


13.10 


11.75 


25.1 


23.9 


21.4 


64.4 






(±0.23) 


(±0.54) 


(±0.65) 


(±1.2) 


(±1.4) 


(±3.4) 


(±4.6) 


a 


[AU] 


0.02225 


0.0641 


0.1286 


0.2699 


0.4929 


1.422 


3.40 






(±0.00035) 


(±0.0010) 


(±0.0020) 


(±0.0042) 


(±0.0078) 


(±0.026) 


(±0.11) 


N 










190 








Span 


[days] 








2428 








rms 


[ms-'] 








1.28 


















1.54 









For that purpose, we added to the minimization an additional 
term corresponding to these proper modes: 

xI^R{u\ + u'^, (9) 

where 7? is a positive constant, that is chosen arbitrarily in order 
to obtain a small value for u \ and U2 simultaneously, without any 
significant decrease of the quality of the fit. 

Using R = 350 we get mj ~ 0.0017 and U2 ~ 0.044 and 
obtain a final - ^ -54, which is nearly identical to the results 
obtained without this additional constraint (R - Q, x^ = 1.50). 
The best-fit solution obtained by this method is listed in TableQ 
We believe that this solution is a more realistic representation of 
the true system than the one listed in Table|4] and we will adopt it 
henceforward for the remaining of the paper Actually, with this 
constraint, the variation of the eccentricity of the two innermost 
planets are now smaller (Fig|8};). 

7.5. Stability of the system 

To analyze the stability of the nominal solution in TableQwe per- 
formed a global frequency analysis (Laskar 1993) in the vicin- 
ity of the six outermost planets (Fig.fTOli. in the same way as 
achieved for the HD 202206 or the HD 45364 planetary systems 
(ICorreia et al.ll2005l l2009h . For each planet, the system is inte- 
grated on a regular 2D mesh of initial conditions, with varying 
semi-major axis and eccentricity, while the other parameters are 
retained at their nominal values. The solution is integrated over 
10 kyr for each initial condition and a stability indicator is com- 
puted to be the variation in the measured mean motion over the 
two consecutive 5 kyr intervals of time. For regular motion, there 
is no significant variation in the mean motion along the trajec- 
tory, while it ca n vary significantly fo r chaotic trajectories (for 
more details see lCouetdic et al.ll2010l) . The result is reported in 



color in Fig.[TOl where "red" represents the strongly chaotic tra- 
jectories, and "dark blue" the extremely stable ones. 

In all plots of Fig (To] the zones of minimal x^ obtained in 
comparing with the observations appear to belong to stable "dark 
blue" areas. This picture thus presents a coherent view of dynam- 
ical analysis and radial velocity measurements, and reinforces 
the confidence that the sub-system formed by the six outermost 
planets given in Table|2]is stable for long timescales. 

Nevertheless, due to the large number of planets in the sys- 
tem, many mean motion resonances can be observed, several of 
them being unstable. None of the planets determined by the solu- 
tion in TableQare in resonance, but some of them lie in between. 
In particular, the pair d and e is close to a 3:1 mean motion res- 
onance, and the pair e and / is close to a 5:2 mean motion reso- 
nance (similar to Jupiter and Saturn). 

7.6. Long-term orbital evolution 

The estimated age of the HD 10180 system is about 4.3 Gyr 
(Table|2]l, indicating that the present planetary system had to sur- 
vive during such a long timescale. We tested the stability of the 
system by performing a numerical integration of the orbits of the 
6 outermost planets in Table|2]over 150 Myr using t he symplec- 
tic integrator SABAC4 of Laskar & Robutel (l200lh with a step 
size of 10"^ years, including general relativity. 

The orbits, displayed in Fig.[TTl remain very stable through- 
out the simulation. We have verified by frequency analysis that 
these orbits evolve in a very regular way. We have also inte- 
grated the full 7-planet system over 1 Myr with a step size of 
10"^ years, without any sign of strong instability, although the 
frequency analysis of the solutions with 7 planets shows that 
these solutions are not as well approximated by quasiperiodic 
series as the 6-planet solutions. This will have to be analyzed 
further as, like in our solar system, the presence of this inner- 




Fig. 10. Global view of the dynamics of the HD 10180 system for variations of the eccentricity and the semi-major axis of the six 
outermost planets. The color scale is the stability index obtained through a frequency analysis of the longitude of the planets over 
two consecutive time intervals. Red areas correspond to high orbital diffusion (instability) and the blue ones to low diffusion (stable 
orbits). Labeled lines give the value ofx^ obtained for each choice of parameters. 



Table 8. Semi-major axes and eccentricity minima and maxima 
observed over 1 Myr in the 7-planet solution of Table|71 
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0.022253 


0.022253 


0.000 


0.082 
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0.064114 


0.064122 


0.010 


0.203 
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0.128536 


0.128626 


0.000 


0.179 
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0.269814 


0.270092 


0.000 


0.156 
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0.492348 


0.493184 


0.023 


0.137 


6 


1.419645 


1.424347 


0.188 


0.242 


7 


3.387207 


3.402716 


0.044 


0.081 



most pl anet seems to be critical for the long-term stability of the 
system (iLaskar & Gastineaull2009l) . 

The fact that we are able to find a stable solution compati- 
ble with the observational data can still be considered as a good 
indicator of the reliability of the determination of the HD 10180 
planetary system. 

Because of the strong gravitational secular interactions be- 
tween the planets, their orbital eccentricities present signifi- 
cant variations, while their semi-major axes are almost constant 
(Table [8]l, which is also the signature that the system is far from 
strong resonances. As the secular frequencies gk (Table|6]l are 
relatively large, the secular variations of the orbital parameters 
are more rapid than in our Solar System, which may allow us 
to detect them directly from observations, and hence access the 



true masses and mutua l inclinations of the planets as it was done 
for the GJ 876 system dCorreia et al.ll2010h . 

7. 7. Additional constraints 

The stability analysis summarized in Fig.[TO]shows a good agree- 
ment between the "dark blue" stable areas and the contour 
curves. We can thus assume that the dynamics of the seven plan- 
ets is not disturbed much by the presence of an additional body 
close-by. 

We then tested the possibility of an additional eighth planet 
in the system by varying the semi-major axis and the eccentricity 
over a wide range, and performing a stability analysis (Fig.fTSll. 
The test was done for a fixed K value {K = 0.78 ms"'), similar 
to planet b. 

From this analysis (Fig.fT2li. one can see that stable or- 
bits are possible beyond 6 AU (outside the outermost planet's 
orbit). More interestingly, stability appears to be also possi- 
ble around 1 AU, which corresponds to orbital periods within 
300 - 350 days, between the orbits of planets / and g, exactly 
at the habitable zone of HD 10180. Among the already known 
planets, this is the only zone where additional planetary mass 
companions can survive. With the current HARPS precision of 
~1 m s"' , we estimate that any object with a minimum mass M > 
10 Me would akeady be visible in the data. Since this does not 
seem to be the case, if we assume that a planet exists in this 
stable zone, it should be at most an Earth-sized object. 
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Fig. 7. Phased RV curves for all signals in the 7-Keplerian 
model. In each case, the contribution of the other 6 signals has 
been subtracted. 



0.6 




2 4 6 8 10 12 14 16 18 20 



time (kyr) 

Fig. 8. Evolution of the eccentricities of planets b (red) and c 
(green) during 20 kyr for three different models. In the top pic- 
ture the initial eccentricity of planet b is set at zero, but mutual 
gravitational perturbations increase its value to 0.4 in less than 
2 kyr (Table|4]i. In the middle figure we included general rela- 
tivity, which calms down the eccentricity variations of the inner- 
most planet, but still did not prevent the eccentricity of planet 
c to reach high values. In the bottom figure we use a model 
where the eccentricities of both planets were previously damped 
by tidal dissipation (Table|7]). This last solution is stable at least 
for 1 Myr. 



We can also try to find constraints on the maximum masses 
of the current seven-planet system if we assume co-planarity of 
the orbits. Indeed, up to now we have been assuming that the 
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Fig. 9. Tidal evolution of the amplitude of the proper modes mi 
(red), M2 (green), M3 (blue), and M4 (pink) resulting from the tidal 
dissipation on planet b with kzlQ- 0.0015. 
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Fig. 11. Long-term evolution of the HD 10180 planetary system 
over 150 Myr starting with the orbital solution from Table|7] The 
panel shows a face-on view of the system, x and y are spatial co- 
ordinates in a frame centered on the star Present orbital solutions 
are traced with solid lines and each dot corresponds to the posi- 
tion of the planet every 10 kyr The semi-major axes (in AU) are 
almost constant, but the eccentricities undergo significant varia- 
tions (Table lU. 




Fig. 12. Possible location of an additional eighth planet in the 
HD 10180 system. The stability of a small-mass particle in the 
system is analyzed, for various semi-major axes and eccentrici- 
ties, and for K = 0.78 ms The stable zones where additional 
planets could be found are the "dark blue" regions. 



inclination of the system to the line-of-sight is 90°, which gives 
minimum values for the planetary masses (Table]?]). 

By decreasing the inclination of the orbital plane of the sys- 
tem, we increase the mass values of all planets and repeat a sta- 
bility analysis of the orbits, as in FigurefTOl As we decrease the 
inclination, the stable "dark-blue" areas become narrower, to a 
point that the minimum of the best fit solution lies outside the 



stable zones. At that point, we conclude that the system cannot 
be stable anymore. It is not straightforward to find a transition 
inclination between the two regimes, but we can infer from our 
plots that stability of the whole system is still possible for an 
inclination of 30°, but becomes impossible for an inclination of 
10°. Therefore, we conclude that the maximum masses of the 
planets are most probably obtained for an inclination around 20°, 
corresponding to a scaling factor of about 3 with respect to min- 
imum masses. 



8. On the properties of low-mass planetary systems 

8. 1 . Dynamical architecture 

The increasing number of multi-planet systems containing at 
least three known planets greatly extends the possibilities to 
study the orbital architectures of extrasolar planetary systems 
and compare them to our Solar System. Although there are al- 
ready 15 systems with at least three planets as of May 2010, one 
should recognize that our knowledge of many of them is still 
highly incomplete due to observational biases. The RV technique 
finds the most massive, close-in planets first in each system, and 
the secure detection of multiple planets requires a large number 
of observations, roughly proportional to the number of planets 
for RV signals well above the noise floor Lower-amplitude sig- 
nals like the ones induced by ice giants and super-Earths require 
a noise floor at or below ~1 m s"' to keep the number of obser- 
vations within a reasonable range. Moreover, the phase of each 
signal must be sufficiently well covered, which requires a large 
enough time span and appropriate sampling. As a consequence, 
the planet detection limits in many of the 15 systems mentioned 
above probably do not reach down to the Neptune-mass range 
yet, preventing us from having a sufficiently complete picture 
of them. Nevertheless, considering the well-observed cases fol- 
lowed at the highest precision, the RV technique shows here its 
ability to study the structure of planetary systems, from gas gi- 
ants to telluric planets. 

The dynamical architecture of planetary systems is likely to 
convey extremely useful information on their origins. The dom- 
inant planet formation scenario presently includes several phys- 
ical processes that occur on similar timescales in protoplane- 
tary disks: formation of cores, preferentially beyond the ice line, 
through accretion of rocky and icy material; runaway gas accre- 
tion on cores having reached a critical mass, rapidly forming gi- 
ant planets; inward migration of cores (type I) and giant planets 
(type II) through angular momentum exchange with the gaseous 
disk; disk evolution and dissipation within a few Myr; dynamical 
interactions between protoplanets leading to eccentricity pump- 
ing, collisions or ejections from the system. It is extremely chal- 
lenging to build models that include all these effects in a con- 
sistent manner, given the complicated physics involved and the 
scarce observational constraints available on the early stages of 
planet formation. In particular, attempts to simultaneously track 
the formation, migration and mutual interactions of several pro- 
toplanets are still in their infancy. Observational results on the 
global architecture of planetary systems may therefore provide 
important clues to determine the relative impact of each process. 

From the observational point of view, the least that can be 
said is that planetary systems display a huge diversity in their 
properties, which after all is not surprising given their complex 
formation processes. Fig.[T3]shows planet semi-major axes on a 
logarithmic scale for the 15 systems known to harbour at least 
three planets. Systems are shown in increasing order of their 
mean planetary mass, while individual planet masses are illus- 
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Table 9. Exponential fits to semi-major axis distributions. 
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Fig. 13. The 15 planetary systems with at least three known plan- 
ets as of May 2010. The numbers give the minimal distance 
between adjacent planets expressed in mutual Hill radii. Planet 
sizes are proportional to log (m sin i). 

trated by varying dot size. The full range of distances covered 
by each planet between periastron and apastron is denoted by 
a black line. The minimal distance between each neighbouring 
pair of planets is also given, expressed in units of the mutual Hill 
radius which is defined as: 
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ai + 02 nil + ni2 



H,M 



3M, 
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(10) 



As shown by [Chambers et al.l (Il996h . the instability 
timescale of a coplanar, low-eccentricity multi-planet system 
is related to the distance between planets, expressed in mu- 
tual Hill radii, in a relatively simple manner, and approximate 
'life expectancies' of planetary systems can be estimated based 
on these separations, the number of planets, and their masses. 
The Chambers et al. simulations do not extend to masses above 
~3 M® or to timescales above 10** years, but a moderate extrap- 
olation of their results shows that for systems with 3-5 plan- 
ets and masses between a few M® and a few Mj, separations 
between adjacent planets should be of at least 7-9 mutual Hill 
radii to ensure stability on a 10-Gyr timescale. These numbers 
should not be taken too exactly since they were obtained assum- 
ing regularly-spaced, equal-mass bodies. They are also not ap- 
plicable to eccentric orbits and dynamical configurations such 
as mean-motion resonances, where stability 'islands' do exist 
at reduced spacings (e.g. GJ 876). However, the global picture 
emerging from Fig. [13] shows that many known planetary sys- 
tems are dynamically ' packed', with little or no space left fo r ad- 
ditional planets (e.g. lBarnes & Ouinnl2004tlBarnes & Ravmondl 
12004 . This result was already noted by several authors, giv- 
ing r ise for example to the 'pa cked planetary system' hypoth- 
esis (iBarnes & Ravmondll200 4'). Here we show that this seems 
to be true also for several low-mass systems, i.e. those which 
do not contain gas giants (or only distant ones), as illustrated 
by HD 40307, GJ 581 and HD 10180. Indeed, several planets 
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23.4 


0.011 
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in these systems are separated from each other by typically less 
than 15 mutual Hill radii. There are still a few 'empty' places, 
however, and further observations will tell if smaller planets are 
hiding between the known ones. 

8.2. Extrasolar Titius-Bode-like laws? 

It is intriguing that many gas giant and low-mass systems seem 
to share the property of being dynamically packed. An attrac- 
tive explanation would be that at each moment of their history, 
many planetary systems are 'saturated' with planets and exhibit 
dynamical configurations whose lifetime is of the same order 
of magnitude as the age of the system. This would point to a 
major role for dynamical interactions in the shaping of plane- 
tary systems, at least since the dissipation of the gaseous disk. 
The observed packing may support the view that close-in low- 
mass systems could be primarily the result of strong interactions 
(collisions and ejections) between several large protoplanets af- 
ter these were brought to the inner regions of the disk through 
type I migration, i.e. after disk dissipation. These systems would 
then naturally evolve towards planets separated from ea ch other 
by a roughly constant number of mutual Hill radii (e.g. iLaskarl 
2000; Ravmond et al. 2009). Since Hill radii are proportional to 
the semi-major axes, the orbital distances of successive plan- 
ets with similar masses will tend to obey an approximate expo- 
nential law, much like the century-long deb ated and polemical 
Titius-Bode law in the Solar System. Indeed. lHaves & Tremaind 
(1998) have shown that any planetary system subject to some 
'radius-exclusion' law such as the Hill criterion is likely to 
have its planets distributed according to a geometric progression. 
iLaskaj yOOO) presents a simplified model of planetary accre- 
tion focusing on the evolution of the angular momentum deficit 
(AMD) of the averaged system. Starting from a given density 
of planetesimals p{a), the final state of the system, defined by 
the end of the collision phase, can be derived analytically and 
the spacing between adjacent planets can be predicted for dif- 
ferent functional forms of p{a). Interestingly, an exponential law 
logfl„ - Cl + C2n is obtained when the initial density p(a) goes 
as a^^^^, while a constant density p(a) yields a semi-major axis 
relation of the form yfa^ - ci + C2n. Studies of the orbital ar- 
rangements in a large number of systems may thus eventually 
allow us to constrain the original density of planetesimals. 

Looking more closely at Fig. [13] we see that exponential 
laws may indeed exist in some planetary systems. However, a 
meaningful test requires that all successive planets have been 
discovered, especially low-mass ones. This is far from being 
guaranteed in the presently-known systems, and we therefore 
limit ourselves to those observed with the HARPS spectrograph, 
and to planets within 1 AU, to minimize observational biases. 
We do not want to speculate on 'missing' planets introducing 
gaps in the Titius-Bode-like relations, since almost anything can 
be fitted to the present, limited datasets if more than two free pa- 
rameters are allowed. As can be seen by eye in Fig. [13] a some- 
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Fig. 14. Fit of exponential laws to semi-major axes as a func- 
tion of planet number for the inner Solar System (black), HD 
40307 (red), GJ 581 (blue), HD 69830 (green) and HD 10180 
(magenta). 



what regular spacing between adjacent planets seems to exist in 
the low-mass systems HD 40307, HD 69830 and HD 10180, but 
less so in GJ 581. Among massive systems, 55 C nc also shows 
a somewhat regular spacing dPoveda & Larall2008l) . but the pres- 
ence of close-in gas giants in this system makes planet-planet 
interactions much stronger and hints at a different formation his- 
tory (e.g. type II instead of type 1 migration) compared to low- 
mass systems. 

Focusing on low-mass systems. Fig. fl4l shows an exponen- 
tial fit a„ - c\ to the observed semi-major axes as a function of 
planet number, starting atn - I. A reasonable fit is obtained for 
HD 40307, HD 69830 and HD 10180, with a relative standard 
deviation of the residuals of 0.57%, 10.2% and 12.0%, respec- 
tively. The fit to the inner Solar System is also shown, with a 
relative standard deviation of 8.0%. The fit to GJ 58 1 is less con- 
vincing, with a dispersion of 21.0%. It is tempting to conjecture 
that there exists an additional body between the third and fourth 
planets in this system, which would make the exponential fit sig- 
nificantly better, and at the same time provide an exciting candi- 
date for a habitable world. Table|9]gives the values of the best-fit 
parameters ci and C2 for each system, together with the average 
planetary masses. Interestingly, a positive correlation between 
C2 and mass may be present, possibly illustrating the fact that 
more massive planetary systems tend to be more widely spaced, 
as would be expected in the context of Hill stability. The ci val- 
ues show how 'special' the inner Solar System is, with the first 
planet (Mercury) very distant from the central star compared to 
the other systems. 

We emphasize that we do not consider these Titius-Bode-like 
'laws' as having any other meaning than a possible signature of 
formation processes. As such, we would expect them to hold 
only in certain types of planetary systems, e.g. close-in, low- 
mass, many-body configurations. The presently-known massive 
systems, on the other hand, likely experienced a more chaotic 
history. Moreover, not all low-mass systems satisfy such expo- 



nential relations (e.g. GJ 581) and the physics of planet forma- 
tion is so diverse and complex that we do not expect any univer- 
sal rule on planet ordering to exist. 

8.3. Formation and evolution 

These emerging patterns, if confirmed by further discoveries of 
planetary systems, may provide clues on how the observed sys- 
tems of close-in super-Earths and Neptunes were formed. These 
systems appear to be quite common, but their formation his- 
tory remains a puzzle. On the one hand, it seems unlikely that 
they formed in situ given the very high inner disk densities 
that would be required. However, little is known about statis- 
tical properties of protoplanetary disks and their density pro- 
files, and this possibility can probably not be completely re- 
jected at this point. On the other hand, such systems may be 
the result of convergent type I migrat ion of planetary cores 
formed at or beyond the ice line (e.g. iTerquem & Papaloizoul 
2007; Kennedv & Kenvon 2008). But how can several proto- 
planets grow to masses in the super-Earth/Neptune range while 
migrating together during the disk lifetime, and end up in a con- 
figuration which is not necessarily close to mean-motion reso- 
nances? Near-commensurability of the orbits would be expected 
according to Terquem & Papaloizou (2007). Loss of commensu- 
rability could occur through orbital decay due to stellar tides, but 
this is probably efficient only for the planets closest to the star 
So this scenario still has difficulties in explaining a system such 
as HD 10180. 

Neve rtheless, as a testable predic tion of this type I migration 
scenario, iKennedv & KenvonI (l2008h suggest that the masses of 
close-in planets will increase as stellar mass increases, and will 
even reach the gas giant range for stars above ~1 Mq. This is 
due the combined effects of increased disk mass, higher isolation 
mass, a more distant snow line and a mass-dependent migration 
timescale to the inner regions that favors more massive planets 
as stellar mass increases. Interestingly, we note that the average 
masses of close-in planets do increase between HD 40307 (M 
= 0.77 Mo, [Fe/H] = -0.31), HD 69830 (M = 0.82 Mq, [Fe/H] 
= -0.06) and HD 10180 (M = I.O6M0, [Fe/H] = 0.08), hinting 
at some correlation between stellar mass, metallicity and masses 
of close-in planetary systems (see also Sect. 18. 41 below). In any 
case, there is great hope that these systems will allow for a much 
better characterization of type I migration in the near future. 

Other mechanisms have been proposed to produce close- 
in low-mass planets, many of which involve the influence 
of (migrating) gas giant(s) furt her ou t in the system (e.g . 
iFogg & Nelsonll2005i IZhou et al.1 120051; iRavmond et al]|2008l) . 
In this context, it is interesting to note that the present RV data 
can exclude the presence of Jupiter-mass objects within ~10 AU 
in the HD 40307, GJ 581, HD 69830 and HD 10180 systems. 
It is therefore unlikely that gas giants played a major role in the 
shaping of these systems. Their absence may actually be the fac- 
tor that favored the formation and survival of many lower-mass 
objects. 

8.4. Correlations witli steilar mass and metaiiicity 

Finally, we may also investigate the impact of stellar mass and 
metallicity on planet formation by further considering the 15 
systems with at least 3 known planets. Fig. [15] shows the total 
planetary mass in these systems as a function of stellar mass 
alone, stellar metallicity alone, and the total amount of heavy 
elements in the star given by MsiaAO^^'^^^\ We note two obvi- 
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Fig. 15. Total planetary system mass as a function of stellar mass (left), stellar metallicity (middle) and the overall metal content of 
the star (right). All 3 quantities are relative to the Sun. The 15 systems with at least 3 known planets are shown. 



ous facts: 1) all very massive systems are found around massive 
and metal-rich stars; 2) the 4 lowest-mass systems are found 
around lower-mass and metal-poor stars. It thus appears that 
both quantities independently impact the mass of formed plan- 
ets. When both effects of stellar mass and metallicity are com- 
bined (right panel), we obtain an even stronger correlation be- 
tween total planetary system mass and total metal content in 
the star. The latter quantity can be seen as a proxy for the total 
amount of heavy elements that was present in the protoplanetary 
disk. These findings confirm previous trends already observed 
for the whole sample of planet-host stars, and are well explained 
by formation scenarios based on the core-accretion model. 



9. Conclusion 

In this paper we have presented a new, very rich planetary system 
with planets ranging from Saturn-mass to Earth-mass, and com- 
prising 5 Neptune-like objects. Long-term radial velocity mon- 
itoring at 1ms ' precision was necessary to detect the low RV 
amplitudes of these planets. The dynamical architecture of this 
system reveals a compact configuration, with planets roughly 
equally spaced on a logarithmic scale, and with significant sec- 
ular interactions. The presence of an Earth-mass body at 0.022 
AU has important implications for the dynamics of the system 
and highhghts the role of tidal dissipation to guarantee stabil- 
ity. Future measurements will allow us to confirm the existence 
of this planet. The HD 10180 system shows the ability of the 
RV technique to study complex multi-planet systems around 
nearby solar-type stars, with detection limits reaching rocky/icy 
objects within habitable zones. Future instruments like VLT- 
ESPRESSO will buUd on the successful HARPS experience and 
carry out a complete census of these low-mass systems in the so- 
lar neighborhood, pushing towards planets of a few Earth masses 
at 1 AU. 

With the advent of the space observatories CoRoT and 
Kepler, low-mass planets have also become accessible to transit 
searches. According to early announcements, the Kepler mission 
will soon confirm what radial velocity surveys are already start- 
ing to find: rocky/icy planets are very common in the Universe. 
The combination of both techniques is hkely to bring rapid 
progress in our understanding of the formation and composition 
of this population. 



The HD 10180 system represents an interesting example of 
the various outcomes of planet formation. No massive gas giant 
was formed, but instead a large number of still relatively massive 
objects survived, and migrated to the inner regions. Building a 
significant sample of such low-mass systems will show what are 
the relative influences of the different physical processes at play 
during planet formation and evolution. 
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Table 1. HARPS radial velocities, CCF FWHM, bisector span and logi?^ measurements for HD 10180. 
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.53884 


0, 


.00055 


7, 


.46672 





.00115 


-0 


.00623 


0, 


.00098 


-4 


.9992 


0, 


.0029 


54733, 


.743640 


35 


.53872 


0, 


.00043 


7, 


.46921 





.00082 


-0 


.00509 


0, 


.00070 


-4 


.9951 


0, 


.0019 


54734, 


,759629 


35, 


.53230 


0, 


,00054 


7, 


.46849 





.00113 


-0 


.00613 


0, 


.00096 


-4, 


.9997 


0, 


.0028 


54736, 


,683562 


35, 


.52815 


0, 


,00057 


7, 


.47328 





.00120 


-0 


.00656 


0, 


,00102 


-5, 


.0024 


0, 


,0032 


54737, 


.736214 


35, 


.53397 


0, 


,00073 


7, 


.47191 





.00162 


-0 


.00589 


0, 


,00138 


-4, 


.9969 


0, 


,0047 


54743, 


.677098 


35 


.54015 


0, 


.00052 


7, 


.47390 





.00108 


-0 


.00517 


0, 


.00092 


-4 


.9993 


0, 


.0027 


54745, 


.663551 


35 


.53661 


0, 


.00047 


7, 


.47070 





.00094 


-0 


.00448 


0, 


.00080 


-5, 


.0042 


0, 


.0023 


54747, 


,822128 


35, 


.52527 


0, 


,00047 


7 


.46858 





.00094 


-0 


.00778 


0, 


,00080 


-4, 


.9985 


0, 


,0023 


54748, 


,786221 


35, 


.52669 


0, 


,00041 


7, 


.46656 





.00078 


-0 


.00511 


0, 


,00066 


-4, 


.9992 


0, 


,0019 


54750, 


,803894 


35, 


.53364 


0, 


,00052 


7, 


.46808 





.00108 


-0 


.00727 


0, 


,00092 


-4, 


.9964 


0, 


,0027 


54758, 


.682427 


35 


.52975 


0, 


.00043 


7, 


.46884 





.00085 


-0 


.00534 


0, 


.00072 


-5 


.0056 


0, 


.0023 


54759, 


.674222 


35, 


.52850 


0, 


.00042 


7, 


.46828 





.00080 


-0 


.00786 


0, 


.00068 


-5, 


.0008 


0, 


.0020 


54760, 


,686923 


35 


.53114 


0, 


,00041 


7, 


.46624 





.00078 


-0 


.00590 


0, 


,00066 


-5, 


.0017 


0, 


,0019 


54761, 


.677044 


35, 


.53778 


0, 


.00043 


7, 


.46567 





.00082 


-0 


.00630 


0, 


.00070 


-5 


.0132 


0, 


.0022 


54763, 


.663210 


35 


.53430 


0, 


.00048 


7 


.46798 





.00097 


-0 


.00800 


0, 


.00082 


-5, 


.0116 


0, 


.0027 


54764, 


,605129 


35, 


.53084 


0, 


,00087 


7, 


.47410 





.00198 


-0 


.00824 


0, 


,00168 


-4, 


.9786 


0, 


,0071 


54765, 


,638628 


35, 


.53012 


0, 


,00069 


7, 


.46935 





.00153 


-0 


.00946 


0, 


,00130 


-4, 


.9948 


0, 


,0049 


54766, 


,631754 


35, 


.53086 


0, 


,00069 


7, 


.46599 





.00153 


-0 


.00433 


0, 


,00130 


-4, 


.9901 


0, 


,0049 


54767, 


.620499 


35, 


.53819 


0, 


.00049 


7, 


.47105 





.00101 


-0 


.00623 


0, 


.00086 


-5 


.0054 


0, 


.0028 


54774, 


.614543 


35 


.54531 


0, 


.00045 


7, 


.47178 





.00089 


-0 


.00691 


0, 


.00076 


-4 


.9899 


0, 


.0021 


54775, 


.656999 


35, 


.54138 


0, 


.00054 


7, 


.47336 





.00113 


-0 


.00571 


0, 


,00096 


-4, 


.9855 


0, 


,0028 


54776, 


,565884 


35, 


.53901 


0, 


,00054 


7, 


.47442 





.00113 


-0 


.00622 


0, 


,00096 


-4, 


.9897 


0, 


,0028 


54777, 


,586673 


35, 


.53941 


0, 


,00057 


7, 


.47079 





.00120 


-0 


.00695 


0, 


,00102 


-4, 


.9872 


0, 


,0030 


54778, 


.583910 


35, 


.54286 


0, 


.00043 


7, 


.47242 





.00082 


-0 


.00615 


0, 


.00070 


-4 


.9890 


0, 


.0019 


54779, 


.561541 


35 


.54268 


0, 


.00050 


7, 


.47296 





.00104 


-0 


.00683 


0, 


.00088 


-4 


.9926 


0, 


.0025 


54780, 


.578359 


35, 


.53982 


0, 


,00047 


7, 


.47084 





.00094 


-0 


.00591 


0, 


.00080 


-4, 


.9912 


0, 


.0022 


54799, 


,586422 


35, 


.52535 


0, 


,00054 


7, 


.47137 





.00113 


-0 


.00551 


0, 


.00096 


-5, 


.0000 


0, 


,0029 


54800, 


,616245 


35, 


.52437 


0, 


,00044 


7, 


.46770 





.00085 


-0 


.00902 


0, 


,00072 


-4, 


.9903 


0, 


,0020 


54801, 


.621649 


35 


.52932 


0, 


.00047 


7, 


.47042 





.00092 


-0 


.00427 


0, 


.00078 


-4 


.9960 


0, 


.0022 


54802, 


.638602 


35 


.53103 


0, 


.00051 


7, 


.47385 





.00106 


-0 


.00727 


0, 


.00090 


-5, 


.0068 


0, 


.0031 
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BID-2,400,000 




RV 




0"RV 


FWHM 


0" FWHM 




BIS 




cbis 


loj 










[k 


ms"'] 


[kms-i] 


[kms-i] 


[kms-'] 


[kms-M 


[kms-'] 


[dex] 


[dex] 


54803, 


.641708 


35, 


.52828 





.00041 


7, 


.46591 


0, 


.00078 


-0 


.00583 


0, 


.00066 


-4 


.9962 


0, 


.0020 


54804, 


,664408 


35, 


.52373 





.00049 


7, 


.46941 


0, 


,00099 


-0 


.00622 


0, 


,00084 


-5, 


.0004 


0, 


,0025 


54805, 


.597303 


35, 


.52220 





.00049 


7, 


.46922 


0, 


.00099 


-0 


.00556 


0, 


.00084 


-4, 


.9989 


0, 


.0024 


54806, 


.581063 


35 


.52765 





.00059 


7, 


.46920 


0, 


.00127 


-0 


.00895 


0, 


.00108 


-5, 


.0034 


0, 


.0038 


54848, 


,523947 


35, 


.53092 





.00050 


7, 


.46837 


0, 


,00101 


-0 


.00581 


0, 


,00086 


-4, 


.9993 


0, 


,0025 


54851, 


,533166 


35, 


.52169 





.00044 


7, 


.46785 


0, 


,00085 


-0 


.00642 


0, 


,00072 


-4, 


.9979 


0, 


,0020 


54853, 


,542410 


35, 


.52842 





.00048 


7, 


.46739 


0, 


,00094 


-0 


.00456 


0, 


,00080 


-5, 


.0000 


0, 


,0023 


55020, 


.896436 


35, 


.54196 





.00044 


7, 


.46884 


0, 


.00087 


-0 


.00573 


0, 


.00074 


-4 


.9911 


0, 


.0020 


55021, 


.904529 


35 


.54027 





.00056 


7, 


.46627 


0, 


.00094 


-0 


.00383 


0, 


.00080 


-4 


.9850 


0, 


.0022 


55022, 


,895816 


35, 


.54093 





.00057 


7, 


.47520 


0, 


.00122 


-0 


.00784 


0, 


.00104 


-4, 


.9874 


0, 


,0029 


55024, 


,901467 


35, 


.53184 





.00050 


7, 


.47442 


0, 


,00104 


-0 


.00552 


0, 


.00088 


-4, 


.9924 


0, 


,0025 


55025, 


,767586 


35, 


.53675 





.00089 


7, 


.47401 


0, 


,00203 


-0 


.00957 


0, 


,00172 


-4, 


.9999 


0, 


,0068 


55037, 


.933241 


35 


.53957 





.00045 


7, 


.46853 


0, 


.00089 


-0 


.00720 


0, 


.00076 


-4 


.9987 


0, 


.0021 


55038, 


.867190 


35 


.53824 





.00044 


7, 


.47193 


0, 


.00087 


-0 


.00558 


0, 


.00074 


-4 


.9991 


0, 


.0021 


55039, 


,919361 


35, 


.53814 





.00196 


7, 


.46885 


0, 


.00459 


-0 


.00664 


0, 


,00390 


-4, 


.9725 


0, 


,0211 


55040, 


,848622 


35, 


.53190 





.00063 


7, 


.47484 


0, 


.00139 


-0 


.00665 


0, 


,00118 


-5, 


.0013 


0, 


,0038 


55041, 


,805672 


35, 


.52934 





.00059 


7, 


.47448 


0, 


.00127 


-0 


.00914 


0, 


.00108 


-4, 


.9981 


0, 


,0034 


55042, 


.876081 


35 


.53238 





.00042 


7, 


.46870 


0, 


.00082 


-0 


.00665 


0, 


.00070 


-4 


.9864 


0, 


.0019 


55044, 


.885876 


35, 


.53455 





.00065 


7, 


.46926 


0, 


.00141 


-0 


.00550 


0, 


.00120 


-5, 


.0119 


0, 


.0039 


55045, 


,847936 


35 


.52777 





.00054 


7, 


.47078 


0, 


,00113 


-0 


.00632 


0, 


,00096 


-4, 


.9952 


0, 


,0028 


55067, 


.829962 


35, 


.53283 





.00050 


7, 


.46967 


0, 


.00104 


-0 


.00583 


0, 


.00088 


-4, 


.9912 


0, 


.0025 


55068, 


.857573 


35 


.53175 





.00059 


7 


.46815 


0, 


.00104 


-0 


.00744 


0, 


.00088 


-4 


.9923 


0, 


.0026 


55069, 


.809858 


35, 


.52925 





.00050 


7, 


.46719 


0, 


.00104 


-0 


.00676 


0, 


.00088 


-4 


.9903 


0, 


.0025 


55070, 


.808393 


35, 


.53095 





.00063 


7, 


.47245 


0, 


.00139 


-0 


.00712 


0, 


.00118 


-4 


.9961 


0, 


.0038 


55071, 


,806534 


35, 


.53285 





.00045 


7, 


.46862 


0, 


,00092 


-0 


.00712 


0, 


,00078 


-4, 


.9956 


0, 


,0023 


55072, 


.870486 


35, 


.53447 





.00044 


7, 


.46680 


0, 


.00089 


-0 


.00722 


0, 


.00076 


-4 


.9980 


0, 


.0024 


55074, 


.843681 


35 


.52982 





.00048 


7, 


.46561 


0, 


.00101 


-0 


.00398 


0, 


.00086 


-4 


.9955 


0, 


.0028 


55076, 


,815394 


35, 


.52740 





.00077 


7, 


.46954 


0, 


,00172 


-0 


.00895 


0, 


,00146 


-4 


.9937 


0, 


,0059 


55095, 


,813607 


35, 


.53214 





.00053 


7, 


.46398 


0, 


,00113 


-0 


.00612 


0, 


,00096 


-5, 


.0037 


0, 


,0030 


55096, 


,775849 


35, 


.52873 





.00052 


7, 


.46821 


0, 


,00108 


-0 


.00639 


0, 


,00092 


-4, 


.9965 


0, 


,0027 


55097, 


.862579 


35, 


.52432 





.00065 


7, 


.47114 


0, 


.00122 


-0 


.00791 


0, 


.00104 


-5 


.0068 


0, 


.0033 


55099, 


.817834 


35 


.52230 





.00047 


7, 


.46958 


0, 


.00094 


-0 


.00540 


0, 


.00080 


-4 


.9967 


0, 


.0023 


55101, 


.809197 


35, 


.53485 





.00093 


7, 


.47209 


0, 


,00212 


-0 


.00525 


0, 


,00180 


-5, 


.0196 


0, 


.0072 


55103, 


.753359 


35, 


.52780 





.00070 


7, 


.47297 


0, 


,00155 


-0 


.00681 


0, 


,00132 


-4, 


.9993 


0, 


,0044 


55106, 


,712448 


35, 


.52879 





.00054 


7, 


.47053 


0, 


,00113 


-0 


.00590 


0, 


,00096 


-4, 


.9983 


0, 


,0029 


55108, 


.743117 


35 


.52738 





.00056 


7, 


.46726 


0, 


.00118 


-0 


.00839 


0, 


.00100 


-4 


.9929 


0, 


.0030 


55110, 


.698889 


35 


.51991 





.00054 


7, 


.46586 


0, 


.00085 


-0 


.00613 


0, 


.00072 


-4 


.9889 


0, 


.0019 


55112, 


.713465 


35, 


.52889 





.00060 


7, 


.46375 


0, 


,00130 


-0 


.00863 


0, 


.00110 


-5, 


.0055 


0, 


,0036 


55113, 


,697280 


35, 


.53144 





.00053 


7, 


.46735 


0, 


,00111 


-0 


.00845 


0, 


,00094 


-4, 


.9932 


0, 


,0028 


55115, 


,751504 


35, 


.52800 





.00068 


7, 


.46998 


0, 


,00151 


-0 


.00836 


0, 


,00128 


-5, 


.0076 


0, 


,0041 


55121, 


.790157 


35 


.53599 





.00057 


7, 


.47652 


0, 


.00122 


-0 


.00850 


0, 


.00104 


-5, 


.0148 


0, 


.0033 


55122, 


.741306 


35, 


.53456 





.00054 


7, 


.47215 


0, 


.00113 


-0 


.00543 


0, 


.00096 


-5, 


.0127 


0, 


.0029 


55123, 


.724617 


35 


.53441 





.00050 


7, 


.47332 


0, 


,00104 


-0 


.00654 


0, 


,00088 


-5, 


.0017 


0, 


,0025 


55125, 


.727881 


35, 


.53788 





.00050 


7, 


.46726 


0, 


.00104 


-0 


.00587 


0, 


.00088 


-5 


.0034 


0, 


.0026 


55128, 


.677184 


35 


.53364 





.00051 


7 


.46875 


0, 


.00106 


-0 


.00650 


0, 


.00090 


-5, 


.0029 


0, 


.0027 


55132, 


,672170 


35, 


.53575 





.00044 


7, 


.46646 


0, 


,00087 


-0 


.00630 


0, 


,00074 


-5, 


.0030 


0, 


,0021 


55133, 


,726967 


35, 


.53641 





.00054 


7, 


.46395 


0, 


,00113 


-0 


.00808 


0, 


,00096 


-5, 


.0145 


0, 


,0030 


55135, 


,624472 


35, 


.54342 





.00045 


7, 


.46779 


0, 


,00089 


-0 


.00853 


0, 


,00076 


-5, 


.0143 


0, 


,0022 


55136, 


.690239 


35, 


.54648 





.00065 


7, 


.46744 


0, 


.00120 


-0 


.00482 


0, 


.00102 


-5 


.0123 


0, 


.0033 


55139, 


.585732 


35 


.53264 





.00060 


7, 


.47179 


0, 


.00130 


-0 


.00494 


0, 


.00110 


-5, 


.0066 


0, 


.0035 


55142, 


,620244 


35, 


.53528 





.00065 


7, 


.47202 


0, 


.00141 


-0 


.00813 


0, 


.00120 


-5, 


.0062 


0, 


,0039 


55143, 


,631659 


35, 


.52923 





.00067 


7, 


.46833 


0, 


,00146 


-0 


.00596 


0, 


.00124 


-5, 


.0010 


0, 


,0040 


55151, 


,694827 


35, 


.52843 





.00048 


7, 


.46810 


0, 


,00099 


-0 


.00656 


0, 


,00084 


-5, 


.0073 


0, 


,0025 


55153, 


.650165 


35, 


.53639 





.00055 


7, 


.46894 


0, 


.00115 


-0 


.00734 


0, 


.00098 


-5 


.0069 


0, 


.0029 


55156, 


.630744 


35 


.52411 





.00106 


7, 


.46449 


0, 


.00243 


-0 


.00498 


0, 


.00206 


-5, 


.0209 


0, 


.0087 


55160, 


,621921 


35, 


.53008 





.00052 


7, 


.46712 


0, 


,00108 


-0 


.00558 


0, 


.00092 


-5, 


.0111 


0, 


,0028 


55161, 


,619721 


35, 


.52798 





.00074 


7, 


.46891 


0, 


,00146 


-0 


.00530 


0, 


,00124 


-5, 


.0035 


0, 


,0041 


55349, 


,912174 


35, 


.53404 





.00062 


7, 


.47053 


0, 


,00137 


-0 


.00604 


0, 


,00116 


-4, 


.9978 


0, 


,0040 


55350, 


.949040 


35 


.52715 





.00070 


7, 


.46860 


0, 


.00155 


-0 


.00416 


0, 


.00132 


-5 


.0109 


0, 


.0048 


55351, 


.941713 


35 


.52050 





.00052 


7, 


.47369 


0, 


.00111 


-0 


.00959 


0, 


.00094 


-5, 


.0109 


0, 


.0031 
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Table 1. Continued. 



BID-2,400,000 




RV 


0"RV 


WHM 


0" FWHM 


BIS 




cbis 


loj 










[k 


tns ] 


ri —It 

[kms 


[kms 


n —It 

[kms 


n —It 

[kms M 


[kms-^ 


[dex] 


[dex] 


55352, 


.935047 


35, 


.52156 


0.00052 


7.47300 


0.00108 


-0.00805 


0. 


.00092 


-4 


.9973 


0.0030 


55353, 


,955486 


35, 


.52173 


0.00057 


7.46834 


0.00125 


-0.00543 


0. 


.00106 


-4, 


.9997 


0.0040 


55354, 


.949932 


35, 


.52597 


0.00061 


7.46902 


0.00132 


-0.00548 


0. 


.00112 


-4 


.9994 


0.0036 


55355, 


.919340 


35 


.52818 


0.00068 


7.47083 


0.00151 


-0.00600 


0. 


.00128 


-4 


.9987 


0.0045 


55357, 


.946100 


3j 


coo C A 


u.uuu /y 


/.469c! / 


A AA 1 nn 
0.001 / / 


A AAA/;'? 

-0.00963 


0. 


.00150 


-5, 


.0072 


A AACT 


55358, 


.928483 


35, 


.52082 


0.00069 


7.47357 


0.00153 


-0.00878 


0. 


.00130 


-5, 


.0043 


0.0047 


55359, 


.943460 


35, 


.52614 


0.00065 


7.47068 


0.00144 


-0.00856 


0. 


.00122 


-5, 


.0108 


0.0045 


55372, 


.929844 


35, 


.53892 


0.00054 


7.46661 


0.00115 


-0.00524 


0. 


.00098 


-4 


.9982 


0.0037 


55373, 
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